Detection of He II A4686 A in 77 Carinae 
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ABSTRACT 

We report the detection of the emission line He II A4686 A in r\ Carinae. The equivalent width 
of this line is ^100 mA along most of the 5.5-yr cycle and jumps to ~900 mA just before phase 
1.0, followed by a brief disappearance. The similarity between the intensity variations of this 
line and of the X-ray light curve is remarkable, suggesting that they are physically connected. 
We show that the number of ionizing photons in the ultraviolet and soft X-rays, expected to 
be emitted in the shock wave from the colliding winds, is of the order of magnitude required to 
produce the He II emission via photoionization. 

The emission is clearly blueshifted when the line is strong. The radial velocity of the line 
is generally -100 Km s _1 , decreases steadily just before the event, and reaches -400 Km s _1 at 
ph= 1.001. At this point, the velocity gradient suddenly changes sign, at the same time that the 
emission intensity drops to nearly zero. Possible scenarios for explaining this emission are briefly 
discussed. The timing of the peak of He II intensity is likely to be associated to the periastron 
and may be a reliable fiduciary mark, important for constraining the orbital parameters. 

Subject headings: stars: early-type-stars: individual (77 Carinae)-binaries:general 



1. Introduction 

Massive stars have strong impact on galactic 
environments. Their evolution, however, is not 
very well known, since their intrinsic parameters 
are difficult to determine. This is even more true 
for the mass, that requires a binary companion to 
be weighted. The Luminous Blue Variable (LBV) 
77 Carinae, believed to be one of the most mas- 
sive stars known from its luminosity, has shown 
signatures of binarity, opening the possibility to 
measure this fundamental parameter. 

The first evidence of binarity emerged from the 
true periodicity in the 5.5 year cycle (Damineli 1996; 
van Genderen et al. 2003; Whitelock et al. 2004; 
Corcoran et al. 2004, in preparation). The sec- 
ond step to unveil the binary nature of the star 
was done when specific models could be calcu- 
lated. Damineli, Conti & Lopes (1997, here- 
after DCL) derived highly eccentric orbits from 
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radial velocities of Paschen lines and suggested 
that X-rays are produced by wind-wind col- 
lision. That particular model had problems, 
since emission lines do not trace precisely the 
orbital motion (Davidson 1997). However, pa- 
rameters similar to that of DCL were used to 
reproduce the main features of the X-ray emis- 
sion in 77 Car by wind-wind collision models 
(Corcoran et al. 2001; Pittard & Corcoran 2002, 
hereafter PC). The orbital parameters, however, 
are still not well constrained. 

The subject of this Letter is to present the un- 
expected detection of the He II A4686 A variable 
emission line in 77 Car and to explore how it is 
related to the binary nature of the star. 

2. Observations and Results 

The data presented here are part of a long-term 
spectroscopic monitoring of 77 Car, started in 1989 
at the Coude focus of the 1.6-m telescope of the 
Pico dos Dias Observatory (LNA/Brazil). Full re- 
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suits, including the analysis of many low and high 
excitation lines, will be reported in Damineli et al. 
(2004, in preparation), from which we are using 
the period length, P=2022.1-d (5.536-yr) and the 
time of phase zero (ph=0.0) of the spectroscopic 
event (JD=2452819.8, 29 June 2003), defined by 
the disappearance of the narrow component in the 
He I 6678 A line. 

Data reduction was done with IRAF package 
in the standard way. The spectral resolution was 
degraded to 0.5A/pixel in order to enhance the 
signal-to-noise ratio. Most of the spectra used 
here have S/N>100 in the stellar continuum and 
some of them twice as much, specially around 
ph=0.2-0.8, when many spectra were co-added. 
We detected an emission line at ~4680-85 A (Fig- 
ure 1), that we suggest to be He II A4686 A. We 
searched for possible transitions in the range 4670- 
4700 A and found none from species that are typ- 
ical of the 77 Carinae spectrum. He II A4686 A has 
never been detected with certainty before, as dis- 
cussed by Hillier and Allen (1992 and references 
therein), who reported an upper limit of lA. The 
line showed up in almost all of our high quality 
spectra, except in a few ones of lower S/N. Para- 
doxically, this line is faint in high excitation phases 
and strong during the event, contrarily to the be- 
havior of the high excitation forbidden lines. In 
Figure 1 we display a sample of spectra collected 
in 2003, ordered by phase of the 2022. 1-d period. 

In order to evaluate the errors, we measured 
several times each line, with different assumptions 
for the continuum. The errors were 20-100 mA in 
the equivalent widths (EWs) and 10-110 Km s _1 
in the radial velocities (RVs) of the line centroid 
(Veen) and are displayed in Figures 2 and 3. Those 
are formal errors; systematic errors may also be 
present, for example introduced by the rectifica- 
tion of the continuum. The fit to the contin- 
uum was done by a 3rd order polynomial, what 
introduces only low frequencies and has little in- 
fluence on scales comparable to the extent of a 
spectral line. Structures in the flatfielding im- 
age could introduce distortions in the continuum, 
although of low frequency, since we used always 
the same observational setup and data reduction 
procedures. The EWs are affected by an addi- 
tional type of systematic error: line blending, that 
prevents the assessment of the local continuum 
and consequently the extension of the line wings. 
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Fig. 1— Sample of He II A4686 A line profiles 
labeled by phase of the 2022. 1-d period. The ver- 
tical dotted line indicates the rest wavelength of 
He II A4686 A. 
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This makes the measurements of EWs and FWHM 
(~500-600 Km s _1 ) somewhat underestimated, as 
is the case for almost all other spectral lines in the 
rich spectrum of r\ Car. The position of the line 
centroid is much more robust and is affected only 
by the S/N ratio, consequently the errors can be 
better judged from the overall scatter in the RV 
curve. 

The slit width was kept fixed at ~1.5" , and this 
seems relevant, since a fraction of the He II emis- 
sion appears to be extended. Our data of 1997 
and 2002, when compared with contemporaneous 
spectra taken with FEROS spectrograph at ESO 
(~3.6" fiber entrance), result in smaller EWs. The 
comparison of our 2003 data with those taken with 
STIS on board of HST (wO.l" slit width) result in 
larger EWs, indicating that the EWs are larger for 
wider slits. This is the opposite of that reported 
by Hillier and Allen (1992) for H I, He I and Fc II 
lines in the Homunculus, which are smaller than 
those in the central object. He II seems to be 
intrinsically in emission near the central source, 
since dust scattering would preserve the equiva- 
lent widths. The line shape and radial velocity are 
unchanged with the slit aperture, indicating that 
most of the emission should arise from the central 
object. In fact, long-slit CCD HST/STIS spec- 
tra (52"x 0.1") collected prior to the 2003.5 event 
(T.R. Gull, private communication) show that the 
He II emission was confined to 0.1" (2-pixel limit) 
in the East-West direction. However, the emission 
could be extended in other directions outside the 
slit. Although this question is relevant, the data 
available to us are insufficient to derive any firm 
conclusion and we will restrict our analysis to the 
homogeneous Brazilian set of data. 

In Figure 2 we see that the strength of the 
He II line is very weak along the 5.5-yr cycle at 
the level of EW«100 mA. It starts strengthen- 
ing around ph«0.8, jumps to EW=873 mA just 
before ph=1.0, and after a brief fading, it rises 
again, reaching a local maximum (EW=300 mA) 
at phwl.04. In Figure 3a we display a zoom in 
the EW curve around ph=1.0, showing a peak 
centered at ph=0.994, a fast drop after ph=1.0 
and a minimum around ph= 1.003. In Figure 3b 
we show the RV curve as measured from the ve- 
locity of the line centroid. The velocity becomes 
more and more negative as the system approaches 
ph=1.0, reaching a sharp minimum (V ce „=— 396 



Kms -1 ) at ph=1.001. After ph=1.002, uncertain- 
ties in RVs increase significantly and we will not 
comment them any further. 

3. Interpretation and Discussion 

The interpretation of the intensity curve of 
He II A4686 A as well as of its radial velocity 
may provide important insights about the na- 
ture and structure of the system. The first thing 
to notice is the strong similarity between the 
He II emission curve and the X-ray light curve 
(Corcoran et al. 2001; Ishibashi et al. 1999). Both 
emissions start to enhance at ph=0.8, peak just 
before ph=1.0, drop to near zero intensity at about 
ph=1.001 and return to normal intensity after 
ph~l.l. A noticeable difference between the two 
light curves is that He II has a much narrower 
and prominent peak (centered at ph=0.994) than 
the X-ray light curve, which has its maximum at 
ph=0.985. The X-ray emission is thought to be 
originated from the shock wave of the colliding 
winds in a binary system. Its "eclipse" has been 
modeled by PC, who derived system parameters 
that will be adopted as a reference frame in the 
following discussion. Given this similarity, it is 
tempting to interpret the He II emission as being 
originated from photoionization by X-ray or by 
the ultraviolet (UV) photons associated to this 
source. 

Near periastron, the X-ray and UV ionizing 
photons are mostly unobserved as they are ab- 
sorbed by the intervening gas along the line of 
sight. However, gas near the shock wave (includ- 
ing the winds of the primary and of the secondary 
stars) may still be directly exposed to this ion- 
ization. The He II A4686 A, at its maximum in- 
tensity, has a luminosity of L# e 77~100 L Q (as- 
suming the somewhat uncertain values of V~7.5 
and Av^6 for the central source). This amounts 
to about 9xl0 46 photons per second. In order 
to produce this emission, a luminous source of 
ionizing soft X-ray and ultraviolet photons is re- 
quired; the peak luminosity of the unabsorbed X- 
rays in the range 2-10 kev is only L;r imo b s =67 L 
(Ishibashi et al. 1999). Detailed hydrodynamical 
numerical calculations by PC show that when the 
velocity of the secondary's wind is large (V2^3000 
Km s -1 ) the energy spectrum behaves like a power 
law, at soft X-rays. From their Figure 3 we derive 
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Fig. 2.— He II A4686 A line EW curve for the last 
three cycles - As are for 1997/98; Ms for 1992; and 
xs for 2003. 
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Fig. 3. — Zoom in He II A4686 A variations around 
ph=0. Upper panel: solid line EW curve; dotted 
lines He II emission models. Lower panel: RV 
curve. Same symbols as in Figure 2 



a spectrum of L £ =2.5xl0 35 e 2,7 ergs s 1 kev 1 
for the system parameters favored by those au- 
thors (M 2 =1CT 5 M yr- 1 and V 2 =3000 Km s" 1 ). 
In the case that this spectrum can be extrapo- 
lated down to 54 cV (He + threshold ionization), 
the total luminosity would be L ion =5200 L Q . This 
extrapolation seems to be reasonable since the en- 
ergy loss by the wind from the secondary star alone 
is i? 2 =7500 L Q . The total number of photons 
emitted by this source is Ni O „=2.5xl0 47 ph s _1 
what amounts to «3 times the observed number 
of He II A4686 A photons. If we add to this the 
energy from the shock of the primary star, this 
number will be even larger. We conclude that near 
periastron the production of ionizing photons is of 
the right order of magnitude to produce the ob- 
served He II emission via photoioniozation. It is 
worth noting that only «1% of the kinetic energy 
from the wind of the secondary star is radiated in 
the 2-10 kev band. Where are the other 99% go- 
ing to? The detection of He II emission may be 
the first hint that they are mostly radiated in the 
UV. 

The fact that the bulk of the He II emission is 
short lived in terms of the EW curve (lasting for 
wl% of the orbital cycle) suggests that its peak in- 
dicates the passage of the periastron - and this is 
quite interesting. This timing is not precisely de- 
termined either by the X-ray light curve (as the 
wind of the primary star is optically thick to X- 
ray absorption) or by the high excitation emission 
lines (that are emitted far away) , although the low 
ionization event is generally believed to be associ- 
ated to the periastron approach. The timing of the 
periastron (ph=0.994) would be the second orbital 
parameter determined with accuracy, besides the 
period. 

A solution to the X-ray light curve was ob- 
tained by PC in which the wind-wind shock before 
eclipse is observed through the opening angle in 
such a way that the X-ray source is seen through 
the non-absorbing wind from the secondary star. 
After the shock front passes through the line of 
sight, the X-rays are absorbed by the wind from 
the primary star. In this way the symmetry of 
the observed emission is broken and one gets a 
post-minimum flux that is lower then the pre- 
minimum one. If we assume that the periastron 
passage occurs at ph=0.994, we are able to fit the 
X-ray light curve with the parameters ecc=0.84 



4 



and the longitude of the periastron of the primary 
star, w=212°. These parameters predict a supe- 
rior conjunction for the secondary star with a true 
anomaly of 0=58° at the phase where our obser- 
vations show a sharp drop in EWs and a gradient 
reversal in the RV curve. Does this suggest an 
eclipse? This question brings us to another one: 
where is the region of He II emission located? 

The large values of the velocity cannot be asso- 
ciated to the shock wave itself as the gas in that 
region has nearly the velocity of the center of mass. 
A similar argument could be made with respect to 
the wind of the primary star. The terminal veloc- 
ity of the primary star is about ^500 Km s^ 1 ; 
at the phase under consideration, the shock wave 
is quite near the primary and the velocity of the 
wind is still far from terminal. In addition, the 
shock wave is symmetric with respect to the or- 
bital plane and the line of sight imposes a projec- 
tion of i~45°. In this scenario, it seems difficult to 
combine the emission and radial velocity behavior 
as displayed in Figure 3. In particular, the radial 
velocity is bluest when the emission has dropped 
to its half maximum intensity, after reaching the 
peak. If we suppose that the emission comes from 
the near side of the primary star, ionized to Hc ++ 
by the secondary or by the shock, one would ex- 
pect maximum blue-shift at maximum emission. 

Perhaps a more promising explanation is that 
the He II emission comes from the wind of the 
secondary star, ionized by the external UV emis- 
sion from the shock wave. This would produce 
an ionization front that divides this wind into a 
zone of He + and a zone of He ++ . The He ++ zone 
would be directly exposed to the UV emission and 
an asymmetric configuration is produced with re- 
spect to the secondary star. This asymmetry will 
cause a blue-shifted emission when the system is 
seen at an appropriate angle. As periastron oc- 
curs at ph=0.994 and the longitude of periastron 
is kj=212°, the ionized He II emitting wind from 
the secondary star should naturally display nega- 
tive velocities when seen from the line of sight of 
the observer. The terminal velocity of the wind 
from the secondary star is about ^3000 Km s _1 
(PC), so one might expect higher velocities and 
line widths than observed. However, the emission 
is likely to come from deep into the wind where 
the density is higher and the velocity smaller, still 
far from the terminal velocity. In addition, the in- 



clination of the system is about i?»45 , which re- 
duces the velocity by about ^cosi. Let us assume 
a highly idealized model of an ionized hemisphere 
with concentric layers of radii r and normal veloc- 
ity law: v=v (l-R2/r) /3 (where R2 is the radius of 
the secondary star). As the emissivity is propor- 
tional to the square of the density, the contribution 
of each layer to the line intensity is proportional to 
^r -2 . The peak of the line profile should, there- 
fore, be emitted at velocities that correspond to 
the innermost ionized zone. We estimate that in 
order to produce the observed RVs (-400 Km s _1 ), 
the innermost ionized zone should occur at r« 4 
R2, for (3=1. At periastron, the distance from the 
secondary star to the innermost zone is 60 R Q , to 
the stagnation shock front is 210 R© and to the 
surface of the primary, 360 R Q . 

The observed line widths (FWHM~500-600 
Km s^ 1 ) require a radius that is somewhat 
smaller. However, as pointed out before, the ob- 
served FWHM is a lower limit because the line 
wings/continuum are difficult to be determined 
with accuracy. If the He II emitting wind is in- 
deed eclipsed, one would expect that when the 
intensity is at its half maximum, half of the wind 
should be eclipsed and the velocity should be most 
negative; this is so because, given the structure of 
the shock wave and its skew angle (see PC), at 
this time the portion of the wind with more posi- 
tive velocity is already eclipsed. This prediction is 
confirmed as both events occur at the same phase 
(ph=1.001). 

The sharpness of the fading phase of the He II 
intensity curve (Figure 3a) is a challenging char- 
acteristic that demands restrictive parameters to 
be modeled. For example, if it is due to an eclipse, 
to be caused by the primary star or by its wind, 
this wind cannot be isotropic with a mass loss 
rate of Mi=2.5x I0~ 4 M yr" 1 (PC); this would 
be optically thick for electron scattering at too 
large a radius. As an alternative, one could sup- 
pose that the Homunculus configuration, with a 
bipolar and equatorial disk could be scaled down 
to the size of the binary system. In this situa- 
tion, the polar wind would be responsible (par- 
tially of totally) for the eclipse of the He II emis- 
sion and the equatorial disk, for the X-ray emit- 
ting shock wave. Such a configuration has also 
been proposed on the basis of other observational 
(Smith et al. 2003; van Boekel et al. 2003) and 
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theoretical (Maeder & Desjacques 2001) consider- 
ations. However, results from spectral modeling 
in the optical (Hillier et al. 2001) and infrared in- 
terferometry (van Boekel et al. 2003) have found 
mass-loss rates of Mi^lO^ 3 M yr _1 . Such de- 
terminations are clearly in contradiction with our 
eclipse hypothesis, unless the mass-loss rate is 
highly anisotropic. 

In the wind-wind collision model, the X-ray 
luminosity is inversely proportional to the separa- 
tion of the stars, Lx a D~ x , so that at periastron, 
X-ray luminosity is maximum (Usov 1992). At 
the same time, the angle subtended by the sec- 
ondary star as seen from each point of the shock 
wave is a D~ x , making the fraction of the UV 
and X-ray photons that are absorbed by the wind 
from the secondary star to go with the square of 
the inverse of the distance between the two stars. 
Therefore one expects that the luminosity of the 
He II line to be inversely proportional to the cube 
of the stellar separation: Lneii ce D~ 3 . Figure 
3 shows that the rising phase and the peak of the 
He II EW curve can be well described by a function 
of the form L HeII a L>~ 3 for ccc=0.82-0.84. This 
agreement, however, does not prove that the He II 
emission comes from the wind of the secondary 
star; any model that provides a dependence as the 
inverse of the cube of the distance would be as 
good as well. 

In this simplified model we suppose that the 
emission is highly concentrated. We should keep 
in mind that there are indications for extended 
emission as well, so that the real picture is actually 
more complicated, but its detailed analysis and 
modeling are beyond the scope of this Letter. 

We thank FAPESP and CNPq for financial sup- 
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